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Abstract 

The small-world networks recently introduced by Watts and Strogatz [Nature 393, 440 
(1998)] has attracted much interests in studying the interesting properties of the networks 
without time-delay. However, a signal or influence travelling on the small-world networks often 
associated with time-delay features which are very common in biological and physical networks. 
We develop an analytical approach as well as numerical simulations to try to characterize the 
effect of time-delay on the properties of small-world networks. An analytical expression of the 
fractal dimension of the small-world networks is given and thus compared with the results from 
numerical simulations. Analysis shows that small-world networks with time-delay generally have 
the multifractals property. 
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1 Introduction 

The properties of complicated networks such as internet servers, power grids, forest fires and porous 
media are mainly determined by the structure of connections between the vertices or occupied sites. 
These networks can be generally represented as a graph network. Obviously, one extreme of such 
a graph is the regular network which has a high degree of local clustering and the average distance 
between the vertices is quite large, while the other extreme is the random network which shows 
negligible local clustering and the average distance is relative small. Recently, Watts and Strogatz 
[1] have presented a model for small- world networks with a high degree of local clustering and a small 
average distance. Such a small- world phenomenon can be obtained by adding randomly only a small 
fraction of the connections, and some common networks such as power grids, film stars networks 
and neural networks behave like small- world networks [2-6]. 

Recently, Moukarzel [3] studied the spreading and shortest paths in system with sparse long- 
range connections by using the small-world model. The spreading of some influence such as a forest 
fire, an infectious disease or a particle in percolating media is studied using a simple rule: at each 
time step, the influence propagates from the infected site to all uninfected sites connected to it via a 
link although this link is not necessary physical. A long-range connection or shortcut can simulate 
the spark that starts a new fire spot, the infect site (say person with fiu) suddenly travels to a new 
place (site), or a portable computer with virus that start to connect to the network a new place. The 
Newman- Watts model [2] and Moukarzel model [3] were mainly focusing the immediate response 
of the network by adding randomly some long-distance shortcuts; there was no time delay in the 
network systems. However, in reality, a spark or an infection can not start a new fire spot or new 
infection immediately, it usually takes some time 5, called ignition time or waiting time, to start 
a new flre or infection. Thus the existing models are no longer able to predict the response in the 
networks or systems with time delay. 

The purpose of the present work is to investigate the small-world networks with time delay 
by using or extending the existing the small-world theoretical framework developed by Watts and 
Strogatz (1998) so as to characterise the fractal dimension of the small-world network with time- 
delay. The present model will generally lead to a delay differential equation, whose solution is usually 
very difficult to obtain if it is not impossible. Thus the numerical simulation becomes essential [1]. 
We will take the analytical analysis as far as possible and compare with the results from numerical 
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simulations. In addition, we shall focus mainly on the 1-D and 2-D networks to study the effect of 
time-delay on fractal dimension and other properties of the small-world networks. 



2 Delay Differential Equation for the Small- World Networks 

Now considering a randomly connected network with n-dimensional lattice [7,8,9,10] and n = 
1,2,3,.... Assuming an influence spreads with a constant velocity u in all directions and a new 
infected spot in the other end of a shortcut will start but with a time delay A. Following the 
method developed by Newman and Watts [2] and Moukarzel [3], the total infected volume V{t) 
comes from two contributions: one is the primary influented volume, and the other contribution is 
the secondary volume by shortcuts. 

The primary volume at time t is the influenced part inside a (hyper) sphere of radiu ut and the 
surface TdiutYn — 1), so the primary volume is TdU^ J* When a primary sphere meets a 

shortcut end (with a density 2p of shortcuts on the network where p is the probability of the longc- 
range shortcuts), a new secondary sphere starts randomly at the other end of the shortcut. This 
secondary volume with a time-delay A is r^M" jJ[2pF(i — C — A)](^^^^^d(. By using a continuum 
approach to the network, the total volume satisfies the time-delay equation of Newman- Watts- 
Moukarzel type 

V{t)=rdu'' f C-'[l+C''V{t-C-A)]dC, n = l,2,3,... (1) 
Jo 

where Td is shape factor of a hypersphere in n-dimensions. The length scale of Newman- Watts type 
[2] can be defined as ^ = l/(2p)^/". However, by slightly rescaling the volume V, a generalized 
definition is more appropriate. That is 

(2pfcn)V"' 

where k being some fixed range. By proper rescaling V by and t by u(^~"r(i(n — 1)!)^/", that 
is 

5 = l/r", r = to(r"rd(n -!)!)!/", 5 = Au(r"rrf(n - 1)!)^/". (3) 
Now we can rewrite (1) as 

Sit) = -^^^^£c-'[l + S{T-C-S)]dC, n = 1,2,3,... (4) 

By diff'erentiating the equation twice, we have the following time-delay equation 

This a delay difi'erential equation, whose explicit solutions are not alway possible depending on the 

initial conditions. For a proper initial condition S(t) = for —6 < r < 0, the solution can be 
explicitly obtained by using the time-forwarding method, 

^(-) = E^^(^' n = l,2,3,... (6) 
Clearly, for r 7^ and ^ ^ 0, the above solution degenerates into 

nj 

given by Moukarzel [3]. For n = 1, S{t) = — 1. For n = 2, S{t) = coshr — 1. However, when 
(5/0, there is no such simple expressions. It can be expected that the time-delay parameter 5 will 
have a strong effect on the evolution behaviour S{t) of an infiuence. 
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3 Fractal Dimension of Delay Small- World Networks 



The evolutionary solution of S{t) with time r and time delay S gives the areas affected by the the 

influence such as fire, infectious diseases and computer virus. One can expect that fractal dimension 
of the small- world networks will give a measure of the network properties as first done by Newman 
and Watts (1999) which does not include the effect of time-delay. Since the velocity of the influence 
travels at a constant speed u, we can rewrite the solution (6) in terms of the distance r = t = At/£, 
where A = u{rd{n — 1)!)^/", we have 

S{r) = Y. J"! ' n = 1,2,3,... (8) 



Since S{r) increases as r increases or S{r) ~ r^, the fractal dimension D can be calculated by 
diflFerential formulae 

z> = M, (9) 

ctlogr 

which can be calculated using (8). We now show how the time-delay affects the fractal dimension 
D of the small- networks. For (5 = 0, we have 



^(-)l^-o = E = ;i E exp[re--/"] - 1, (10) 

j=i j=o 



so that fractal dimension becomes 

^. r{lE;Co^exp[^+re^^-^V"] 
^E;=d^xp[^+re'2-i/"]-l 

For n = 1, this becomes D = re^/{e^ — 1). Clearly, when r — >■ 0, the limit of D leads to one. One 
the other hand, when r — >■ oo, the r— term dominates, so we have D ^ r = At/^. For 5^0, there is 
no general explicit asymptotic expressions and numerical evaluations can be easily done. For small 
time delay and shorter length scale (r = At/S, <C 1), the fractal dimension remains nearly a constant. 
For r 1, the fractal dimension increases essentially linearly with the radiu r, thus the dimension 
of a small- world netword depends on the length scale on which one looks at it. The dependence of 
the fractal dimension D on the lengthscale ^ through r (since r and r is rescaled by the length scale 
^) and the time-delay S suggests the multifractal features of the small- networks. 

On the other hand, for a network of finite size, we can change the shortcut probability p to 
modulate the effective connectivity and fractal dimension. For a very small p 0, the effective 
length scale ^ — >■ oo, so any finite-sized network satifies the condition r ^ 1, and thus the fractal 
dimension D — > 1. As p increases, the length scale ^ reduces, so that the fractal dimension increases 
and the network is closely interconnected with increasing efficiency. 

In order to compare the analytical solution with numerical simulations, we use the simulation 
method given by Watts and Strogatz (1998) and Newman and Watts (1999). The numerical simu- 
lations for the one-dimensional n = 1 and the comparison of fractal dimension with the analytical 
solution (8) is shown in Figure 1 where the network size N = 500, 000, ^ = 500 (or p = 0.002 and 
k = 2) and 5 = 0, 1,5, 10. We can see that numerical results (marked with o, o etc) arc in good 
agreement with the analytical solution (solid lines) by using (8). As expected, for r ^ 1, both nu- 
merical results and analytical fractal dimension approach n = 1. The time-delay has a very strong 
effect on the fractal dimension D of the small- world networks. For S ^ 1, D is substantially reduced 
compared with the one without time-delay 6 = 0, and D is quite near n for most of the region. 
The large time-delay actually makes the network become a larger world network, and the inffuence 
slowly spread in more localized areas. 

The transition from a very small to higher value of p or very long to moderate length scale 
^* = l/{2pkny^"' suggests that the increasing shortcut probability p will greatly increase the fractal 
dimension D of the small-world network. This has very important implications for real world net- 
works. For a small-world network such as the internet, road networks, and business partnerships. 
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the introduction of a small fraction of long-range shortcuts or connections, will essentially make the 
network behave more effectively since the effective scale is becoming smaller. Real networks have 
always some time delay that usually makes the network larger than networks without time delay. 
Thus there is a trade-off between these two competitive factors and the proper network design is 
very important to ensure the efficiency of real networks such as the internet as well as business 
partnerships. 

In summary, numerical simulations and analytical analysis for small-world networks with time- 
delay show that the time-delay parameter S has a very strong effect on the fractal dimension and 
other properties such as the speed and saturation time of the delay networks. A small-world network 
can become larger if sufficient time-delay in the system response is introduced. On the other hand, 
in order to make a large-world network transfer into a small- world, a slightly higher probability p of 
long-range random shortcuts are necessary compared with the one without time-delay because this 
essentially reduces the characteristic length scale ^ = 0{l/p) quite significantly. Further studies of 
the dynamics of the delay small- world networks are now in progress. 

Acknowledgement: The author would like to thank the referee (s) for their very helpful com- 
ments which have greatly improved the manuscript. 
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Figure 1: Fractal dimension of small-world networks with time-delay 5 = 0, 1,5, 10 for a network 
size = 500, 000, ^ = 500 (or p = 0.002) and n = 1. Numerical results (marked with o, o etc) agree 
well with analytical express (solid) calculated from (9). 
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